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O^ ' OfF-resonance conductance through weakly coupled quantum dots ("valley conductance") is 

^^ ' governed by cotunneling processes in which a large number of dot states participate. Virtually the 

same states participate in the transport at consecutive valleys, which leads to significant valley- 
valley conductance correlations. These correlations are calculated within the constant interaction 
model. Comparison with experiment shows that these correlations are less robust in reality. Among 
the possible reasons for this is the breakdown of the constant interaction model, accompanied by 
OO ' "scrambling" of the dot as the particle number is varied. 

(N 
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I. INTRODUCTION 



One of the interesting aspects of the physics of quantum dots is the mechanism of cotunneling ^ which governs 
transport through the quantum dot away from resonance (" conductance valley" ) . Such a mechanism, which usually 
"ti ' gives a small (as compared with resonance values) yet significant contributions to the conductance, consists of the 

a use of a large number of virtual dot states, which, due to high electrostatic energy, are classically forbidden. As one 

varies an external parameter (e.g. applied magnetic field), these are virtually the same (albeit possibly modified) 
'"^ ' states that contribute to the transmission, giving rise to significant conductance correlations H. In a recent work 0] 
^ ' we have pointed out that such correlations show up in the transmission phase as well. Here we study the conductance 
O ' at different valleys (corresponding to a different number of electrons on the dot, a parameter which is controlled by 
I ^1 - an applied gate voltage) and calculate the ensuing conductance correlations. Our formalism, outlined in the present 
and throughout the next section, follows that of Aleiner and Glazman B]. We evaluate the significant valley-valley 
CN ■ conductance correlations within the constant interaction model, as function of the strength of the interaction, location 
^ ■ in the valley and temperature (Section III and Appendix A). A simplified toy model which, in our opinion, captures 
r^ much of the pertinent physics, is presented and studied in Appendix B. Comparison with the experiment (Section IV) 
reveals that in reality these correlations are less robust. Possible reasons for this are discussed in Section V, where 
we stress the likelihood of the breakdown of the constant interaction model, accompanied by " scrambling" of the dot 
states as the particle number is varied. 
^\ ' We consider transport through a quantum dot weakly coupled to reservoirs by high tunneling barriers MM- The 

On • chemical potential fi of the quantum dot can be tuned by an applied gate voltage Vg. Because of the small size of 
■4-^ ' those structures charging effects are important; the weak coupling to the external reservoirs (leads) implies that, in 
3 . general, there is a definite number of electrons on the dot. The conductance of the quantum dot as function of Vg 
^ ■ shows nearly equally spaced peaks. These resonances are due to sequential tunneling of electrons through the dot. 
', ' At these values of Vg the ground state energies of the system with the dot containing N and iV + 1 electrons are 
^ • degenerate and the occupation number can fluctuate. In the "valleys" between the resonances the conductance is 
Q [ strongly suppressed due to the charging energy: it costs energy to add or remove electrons to or from the dot, and 
CJ . the number of electrons on the dot is fixed to an integer value. 

There is however a residual conductance in the valley between the peaks which is caused by a different transfer 
process. An electron or hole tunnels virtually from one reservoir to the other through energetically forbidden states. 
p\^ ' Since the particle tunnels coherently through the whole structure this transfer mechanism is known as cotunneling 
Jh I W- In the following we will focus on this regime away from the resonances. 

We model the quantum dot coupled to the reservoirs by the usual Hamiltonian 

H = H^ + H" + H'^° + H'^ , (1) 

H'^-^ =j:e,at^-Kt^^\ (2) 

k 



H^ = Y. ^^4«fe + ^•'^- + L^R, (3) 



i 

jjL,R (jggQj-ibe the reservoirs to the left and right of the quantum dot, H^ represents the tunneUng of electrons in and 
out of the quantum dot, and H^^ describes the states of the isolated quantum dot including the electron-electron 
interaction. We use the simplifying "constant interaction" model which asserts that the total energy due to the 
electron-electron interaction solely depends on the total number of electrons on the dot (N is the number operator). 
The charging energy is given by U — e^/C invoking the electrostatic energy of a classical capacitance C. The coupling 
strength of level j of the quantum dot to the leads is characterized by the tunnehng rates Tj = 2n/h J2k l^j.k P'5(i?— e^). 
For the calculation of transport properties we need the retarded Green function of the quantum dot coupled to the 
leads, 

Gl':;{t) = {[c,{t),c]{0)U). (5) 

In the regime of weak coupling F ^ kT, A (temperature T, mean level spacing A) we approximate G"'®* by 

(»j)jv . 1 - (nj-)jv 



gj:^ = E ^ 



NOjk 



N=0 



E-{ej- fi + U-{N- 1)) + iVj E-{ej- ii + U ■N) + iV^ 



(6) 



with {..)n — tiN exp~(3H^^ .../tr^ exp~f3H'^^ the thermal average with N electrons. Here A is the single particle 
level spacing. The probability to find N electrons on the dot is given by Pn — {^)n /J2m(^)m- Away from the 
resonances eqn. (|g) describes the elastic cotunnling, i.e. the virtual tunneling via a single state j of the quantum dot. 
This is the dominant process in the regime kT < y/UA [Q which we are focussing at. In addition, (^ also describes 
the dynamics in the vicinity of the resonances on the same level as an master equation approach of ref . |q] . 

For the Green function G"°' we note that in the cotunneling regime away from the resonances there is an integer 
number N of electrons on the dot, and Pjv = ^n n ^^ (h)- Hereafter we attach an index TV to a valley, corresponding 
to the number of electrons on the dot over that range of Vg . The resonance separating the valleys iV — 1 and N will 
be denoted by (TV — 1, TV). It is convenient to shift the reference point of the chemical potential /x to the resonance 
{N — 1, TV) {n —^ fi — eN — U{N — 1)), and make use of the parametrization of n by x: for the valley TV the parameter 
a; ^ corresponds to the right of resonance (TV — 1,TV), i.e. the point where the energies of the system with TV — 1 
electrons and TV electrons on the dot are practically degenerate; a; ^- 1 corresponds to the left of resonance (TV, TV + 1). 
A complete description of ^ is given by the two variables (TV, x) . 

Eqn. (g) contains canonical occupation numbers. In the following we will use the grand-canonical occupation 
numbers, i.e. the Fermi functions, instead. We account for the difference by using an effective inverse temperature 
Pgc ~ /5(1 + /5A/4) which is justified for /3A < 1. We then obtain with ui = e^ 

EAT — W + XU UI — £N + (l — X)U 

Since we are in the cotunneling regime (large denominators), we neglect the tunneling rates iV in the denominators 
in eqn. m). The first (second) term describes the occupied (empty) states with energies smaller (greater) than e^r. 
The minimum energy necessary for hole (particle) transfer is -E^ = xU {Ep = (1 — x)U) with _E^ + Ep = U . We 
should keep in mind that eqn. (uu is valid for x sufficiently away from the resonances at a; = 0, 1. Since the number 
of electrons TV starts to fluctuate on an energy scale kT around the resonances the range of validity of eqn. (0) is 
kT/U <x<l- kT/U. 

II. FORMALISM 

We will briefly review the formalism of Aleiner and Glazman |2] for the calculation of transport properties of a 
disordered quantum dot in the cotunneling regime. Hereafter "disorder" should be understood as caused either by 
impurity scatterers or by shape irregularities in the confining potential. In both cases the statistical properties of 
the single-particle levels are described by random matrix theory (RMT) (for a review see Iq]): over energy intervals 
smaller than the Thouless energy Eth, the tunneling matrix elements and the level spacing are strongly fluctuating. 

We start with the expression for the transmission amplitude t{E) containing the retarded Green function of the dot 
G and the tunneling matrix elements V, ' , 



^(^)=E^'^(^)^''*(^)Gf(^)' (8) 



where Vj" = Vh,^-, y/2TTp^{E) with p^ the density of states in the left reservoir. Since in the elastic cotunneling 
regime the electrons or holes are transferred through one single dot state, only the diagonal elements in (^ have to 
be accounted for. The tunneling matrix elements l^fc(_B) in the Hamiltonian (^ contain the overlap of the lead wave 
function with the dot wave function at the barrier |9| . Wc assume that the tunneling matrix elements factorize into a 
part describing the tunnel barrier with its penetration factor, and a dot part. The dot part involves the wave function 
of state j at the barrier, ipj{R). Thus, we have for the left barrier V^^,,^-. = V^i/jj{R^). The separation is justified 
when the properties of the barrier only weakly change on the typical energy scale, U, of the system. Eqn. (1q) then 
becomes 

t{E) = 27Ty/J/^V^V'^*Y,^j(R^WAR'^)GfiE) 

j 

= 2n./^^V'^V^* /dc^ V,(i?^)^*(i?«)<5(u; - e,)GT{E) . (9) 

with p^'^ the density of states in the leads. For the last step in (pi) it is essential that G"^^ only depends on j via 
ej , cf . (0) . Aleiner and Glazman identified the sum over j as the local density of states of single-particle levels and 
expressed it in terms of Green functions of the non-interacting dot 



G^{R^,R")-G^{R\R'')=J2 



V',(i?^)^*(i?«) Vj(i?^)V'*(i?^) 



uj — ej — iS uj ~ ej + iS 

= 2«^^,(i?^)^;(i?^)<5(c^ - £,). (10) 



Impurities on the quantum dot only affect the single-particle states. With (hO|) the transmission amplitude is expressed 
as a convolution of two types of Green functions, 

tiE) = 2^v/77^F^F^* J ^^{Gt{R\ i?^) - GS{R'^,R"))GT{E) . (11) 

One term in the integral is the cotunneling Green function G"'''* of a quantum dot with interaction (but no explicit 
dependence on disorder dependent quantities); the other term includes the advanced and retarded Green functions 
G^ and G^ of the system which account for disorder, but do not include interaction. We would like to stress that this 
separation into interaction and disorder dependent terms respectively is only possible due to the constant interaction 
model. 

The transmission amplitude t(E) is linked with the linear conductance Q through the quantum dot by the Landauer 
formula [noljO]. Off the resonances Q = e"^ / h\t{EF)\'^ even for finite temperature. Averaging over disorder only 
affects the single-particle Green functions G^'^ defined for non-interacting electrons. For weak disorder, kpl ^ 1, a 
diagrammatic expansion can be established |jl3]. The averages (G^G^) and [G^G^) reduce to the product of the 
average of the single particle Green functions; as long as the dot's size L — \R^ — R^\ is larger than the mean free 
path I such terms may be neglected. The remaining average {G^G ) can be expressed in terms of Diffusons and 
Cooperons, 

(G5„B,(i?^i?^)G4,B.(^^,i?'')> = 27rp^<^l^J„ (12) 

(G« ,B,(i?^i?'')G^.,s.(^^^^)> = '^^P°C!^IL% (13) 

where p^ is the density of states on the dot. Diffusons (Cooperons) depend only on the difference (sum) of magnetic 
fields Bi,B2, and are the solution of diffusion equations. In momentum space {R^ <-^ fc, R^ ^^ k) these solutions read 



where q = k — k and Q — k + k. The area of the quantum dot is denoted by S and r2+ = B'^EtkS^ /^^ with $o = hc/e 
the flux quantum. The Cooperon contribution is maximal for (5 = 0, i.e. for the coherent backscattering process 
k = —k. In contrast to the Diffuson, the Cooperon is suppressed for a magnetic field greater than a characteristic 
field Be ~ ^o/S. The lowest non-vanishing eigenvalues of the diffusion operators are of the order of the Thouless 
energy Eth- For energies E < Eth only the zero- frequency mode must be retained which corresponds to spatially 
homogeneous solutions in ( J14| ) and (jl^) {q,Q -^ 0). 

We focus on the regime A <^ Eh,Ep < Eth- The first inequality assures that we are allowed to use the diffusion 
approximation. Since the relevant energy scale (charging energy) is smaller than the Thouless energy Eth which is 
the case for a dot with not too strong disorder or boundary deformations, we only have to consider the spatially 
homogeneous zero mode q -^ Q. The sum of Diffusons in momentum space then becomes 

S^^ 5"-! 

V^{q)+V^^[-q) = 



—iuj + Dq^ ito + Dq^ 
= -iS-^ { ^ ^) "^ 2^S5{uj). (16) 

Aleiner and Glazman have obtained the full distribution function of Cy [p| and have found for the variance 

r /r2\ /r\2 / 2(e)' for 5 = ,^_. 

varg = {g)-{g) = | ^^^^ forB»i?, • (1^) 



III. COTUNNELING CONDUCTANCE CORRELATIONS 

In the conductance-conductance correlation for different magnetic fields but the same chemical potential was 
studied. Here, we will calculate {g{Ni,xi)g{N2,X2)) between different values of the chemical potential (A^i,a;i) and 
(A^2, X2) but fixed magnetic field for the two cases B = and B ^ Be. We are especially interested in the correlation 
function for N2 = Ni + AN and xi = a;2_, i.e. for different valleys but at the same position within the valleys. In 
order to make contact to the experiment |l3] we consider the following normalized autocorrelation function 

r.,.r ., , {Q{Nu^i)Q{N2,X2)) ~ {g{N^,x^)){g{N2,X2)) .^^. 

C{Ni,xuN2,X2)^ , , (18) 

^vary(A'^i, a;i)^vary(iV2, 2:2) 

for U < Erh- We expect to find strong correlations since from one valley to the next only one level changes its 
character from occupied to empty (or vice versa) Q. For instance, the first empty state in the A''*'' valley is the last 
occupied state in the (A^ -|- 1)*** valley. 

We start as shown in the previous section by separating interaction from disorder 

{g{Ni, xi)g{N2, X2)) = K.^ dijJidijj2dijJzduJi'K 

G-*(A^i, xi)GLt(A^i, a;i)*GLt(A^2, X2)G-'(A^2, 0:2)* x 

( {Gt, (i?^ i?^) - G^^ (i?^ i?^)) (G^, (i?^, i?^) - Gt, (i?^, i?^)) X 

(G:^3(i?^i?^)-G^3(i?^i?«)) (G5,(i?«,i?^)-G:^,(i?^,i?^))) . (19) 

with /C = e^ /h\V'"\'^\V^\'^ p^ p^. Since the energy scale {U) involved in this problem is much larger than the mean 
level spacing (A) the average of the product of four Green functions can approximately be decoupled into a product 
of pairwise averaged Green functions y] . In a diagrammatic expansion this approximation corresponds to keeping a 
certain class of diagrams which also turn out to be relevant for the universal conductance fluctuations il4,pql . As in 
the previous section we will neglect terms like {G^G^) and {G^G^). All possible pairings of the remaining averages 
in terms of Diffusons and Cooperons are given by 

((G^ - G")iG^ - G"'){G'^ - G"){G'^ - G^)) = 

(I?(j2 _i^i (i? ,R ) +'D^-^^^^{R ,R ■)) {V^^^^.^{R ,R ) +'D^^^^^{R ,R )) + 

{'^ui-L.jiiR ,R ) +T^uji-uji{R i^ )) (^W2-W3(-R iR ) +'D^^^^^{R ,R )) + 

[Cuj3-ui{R ,R O+Ci^i-wal-R ,R )) \Cuj2-ui4.{R i^ ) +C^^^^^{R ,R )) . (20) 



Due to the Cooperons, the conductance-conductance correlation function (19) depends on the magnetic field. Using 
(P0|), eqn. ([l9| ) can be expressed as 

{g{Nuxi)g{N2,x2)) = {g{Nuxi)){giN2,x2))+ 

2^p^/C J dujiduj2GT,{Ni,xi)G''f^{N2,X2r {V^,-u., {R"", R"") + V^,-^, {R'',R")) 
2^p^/C j duJiduj2GT,{N^,xi)GT,{N2.x2) (C,-^,(i?^,i?«) +C,-^,(i?^,i?^)) 



(21) 

For a strong magnetic field the last term involving the Cooperons vanishes. For U < Exh we only have to consider 
the spatially constant zero mode of the diffusion operator, and the sum of the two Cooperons also tends to 5{uji — UJ2) 
for B = 0, as in (|lj). Thus the conductance correlation function becomes 

{g{N^,xi)g{N2,X2)) ^ {g{N^,xi)){g{N2,X2)) + a{^n''p'']CS-^f\A{NuXi,N2,X2)W 



where A{Ni,xi,N2tX2) 



dLuGTiNi,x,)GT{N2,X2) 



(22) 



The parameter a depends on whether time reversal symmetry is broken or not: a = 2 at -B = 0, while a = 1 for 
B ^ Bf.- Thus the conductance correlation function (18) becomes 



GiNi,xi,N2,X2) 



\AiNi,xi,N2,X2)\^' 



\A{Ni,xi,Ni,xi)\\A{N2,X2,N2,X2)\' 



(23) 



We note that G is independent of a, i.e. possesses a universal form for both conserved and broken time reversal 
symmetry. The integral for A at finite temperature can be evaluated and the result is given in Appendix |A|. At low 
temperature A/3 3> 1 the integration yields with d ~ e^r^ — e^q^ positive 



A(iVi,a;i,iV2,a;2) 



1 



d+ (.T2 - X])U 



log 



d^X2U 
xiU 



1 



d+ {x2 — xi)U 



log 



d+{l-xi)U 

{l-X2)U 



1 



log 



d + X2U 



log 



d+(l -xi) 



d + U +{x2-xi)U \ ° X2U ° {l-xi)U 

The autocorrelation function ([l8|) in this limit finally becomes 



(24) 



G{Ni,xi,N2,X2) = .Tia;2(l -a::i)(l - X2) 

1 



1 



d/U + X2 — X- 



■log 



{d/U + X2){d/U + l-xi) 



a;i(l - X2) 



■log 



{d/U + I - x,)id/U + X2) 



(25) 



d/U + 1 + X2 — xi ° (1 — a;i)x2 

The energy difference d between the single-particle levels N2 and A^i can be expressed in terms of the mean level 
spacing A, d/U = {N2 — Ni)A/U — {N2 — Ni)S where S — A/U. As a special case (and for the sake of simplifying 
the algebraic expressions), we set cc = xi = 0:2 and obtain for the conductance autocorrelation function G for different 
valleys 

G{Ni,x,N2,x) =x^{l-xf X 



1 



■ log(l 



{N2 - N,)5 



)(1 



{N2 - Ni)S 



(26) 



_{1 + iN2 - Ni)6)iN2 - Ni)6 '"^'^ ' X ''' ' 1 

Fig. |l|(left) shows G as function of A^2 — -^1 for U/A ~ 40 and different values of x. As expected the correlations 
decay slowly with valley number. Away from x — 1/2, G falls off more rapidly because the nearest occupied/empty 
level contributes more strongly to the conductance. Fluctuations in this level from one valley to the next reduce the 
correlations. Fig. |l| (center) demonstrates that the scale for the decay is set by U/A. The slope of C at A^2 — A^i =0 
at X = 0.5 is given by —AA/U. 

Next we turn to the correlations between Xi and X2 in the same valley A^i = A^2- The autocorrelation function ( |25| ) 
then becomes 

2:1X2(1 -xi)(l -.X2) ^ ' ^ 



G{N,xi,N,X2) 



log 



X2 1 - Xl 



(27) 



(X2 — Xi)2 \^ ° Xl 1 — X2, 

We note that G is independent of the ratio U/A. Fig. |] (right) shows G vs. xi and X2 fixed. At xi — X2 the 
autocorrelation function becomes maximal, and it falls rapidly off for xi — > 0, 1. We should bear in mind, however, 
that (|^) breaks down near the resonances, since it neglects fluctuations in the particle number of the dot. 




FIG. 1. Conductance autocorrelation function C for a;i = 2:2 = a; vs. valley number A'^2 — Ni. Left: For U/A = 40 and 
X = 0.5 (solid), 0.2 (dotted) and 0.1 (dashed). Center: For x = 0.5 and U/A = 40 (solid), 30 (dotted) and 20 (dashed). Right: 
C for A'^i = N2 vs. relative position in the valley xi and X2 fixed at 0.5 (solid), 0.7 (dotted) and 0.9 (dashed). Note that ( [2 7] ) 
is not valid in the vicinity of order kT/U around xi = and a;i = 1. 

IV. COMPARISON WITH AN EXPERIMENT 



Measurements of the inter- valley conductance of semiconductor quantum dots were first reported in u^. The 
experimental set-up was similar to other experiments reported in ref. ||lq,0| . The dots were all in the ballistic regime 
and an irregular shape of the confining potential renders the motion of the electrons chaotic, so that RMT describes 
the statistics of the mesoscopic fluctuations in the sample. Additional gate electrodes could distort the shape of the 
dot, and allow an ensemble average to be obtained from a single sample. The experiment focused on the conductance- 
conductance correlations in the elastic cotunncling regime between different magnetic fields but at the same chemical 
potential, as calculated in |^. Since kT < ^/UA the observed cotunneling was indeed elastic. In ref. [|3[ the dot was 
more strongly coupled to the leads (barrier conductance Q^'^ ^ e^//i) than in ||l6| , [l7| , since otherwise the cotunneling 
current would be very small and difficult to measure. Thus the tunneling rate T ^ 0.7A and the averaged conductance 
at the peaks is eight times larger than the cotunneling conductance in the valleys {{Q) ^ 0.05e^/h). 

Fig. shows a comparison of C vs. N2 — Ni between theory and experiment (ref. [13[). The quantum dot on the 
sample under study had the parameters kT ^^ 9^eV, A = Ib.SfieV and U — 410/ieT^. Experimental data for the 
autocorrelation function C (Iq) in the middle of the valleys a;i = a;2 = 1/2 are only available for three neighboring 
valleys and for U/A — 25.9. The lines are our results for three different temperatures. We observe that in the 
experiment the correlations are more strongly suppressed than predicted in theory. 



V. DISCUSSION 



There are two main messages arising from the present analysis. On the theoretical side, our analysis underlines 
the significant correlations of the transmission (i.e. conductance) through a quantum dot as an external parameter is 
varied. In our case the quantum dot is weakly coupled to the leads (external reservoirs) and the external parameter 
which is varied is the gate voltage, affecting the number of electrons which reside on the dot. All this results in 
valley-valley correlations of the conductance in the cotunneling (far-from-resonance) regime. Previously, correlations 
have been found in the original analysis of Aleiner and Glazman (intra-valley, as function of the magnetic field) , in 
ref. |l8[ (valley-valley) considering the differential capacitance of the dot instead of the cotunneling conductance, and 
in a recent analysis on transmission phases |3|. The mere existence of correlations as found in the present study is, 
therefore, not totally surprising, although the details (dependence on the location within the valley, on the interaction 
strength-charging energy, temperature etc.) are certainly different. 

The second message contains, in our opinion, a much more intriguing element. The correlations found in the 
experiment seem to be less robust than our theoretical expressions suggest. There might be several possible reasons 
for this. Some are obvious while others are more subtle and may give rise to some intriguing physics. 

(i) The precise value of the electron-gas temperature in the experiment. In M] the temperature was estimated to 
be kT ~ 9^eV which corresponds to /3A — 1.4. Modest deviations from this value are not going to produce 
agreement with theory. 
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FIG. 2. Comparison between experiment (filled circles) and theory for the autocorrelation function C at x = 1/2 for several 
neighboring valleys {U/A — 25.9). Solid line: obtained from (GGf) at /3A S> 1; dashed line: obtained from the expressions for A 
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from appendix ^at Pec A. = 4 which corresponds to /3A — 2.47; dotted: /3gcA = 1.89 corresponding to /3A — 1.4. 



(ii) Approximating the canonical distribution function by a Fcrnii function with an effective shifted temperature 
(cf. the comment that precedes eqn. (0), see also ref. 0) is asymptotically justified for kT ^ A, which, for 
/3A = 1.4 is not quite the case. However, given the other curves in Fig. 0, a more accurate treatment of this 
point is not going to cure the problem. 

(iii) In our derivation of the correlation function C we have considered the case U < Eth ■ From the technical point 
of view, this allowed us to retain the "zero-mode" of the Diffuson propagator only. In the experiment, the 
Thouless energy was estimated to be Exh = ISOfxeV, which is smaller than the charging energy oi U — 410/iey. 
Qualitatively one can expect that the inclusion of non-zero modes would imply deviation from the constant 
interaction model. In that case the addition (removal) of an electron is likely to change the effective (interaction 
induced) potential landscape felt by the other electrons, facilitating more efficient "scrambling" of the dot (see 
below), and suppressing the correlations found here. However a quantitative statement in this regard calls for 
a detailed analysis, not included in the present work. 

(iv) The Hamiltonian of an interacting electron system (such as a quantum dot) includes terms other than the 
constant interaction U . Varying the gate voltage is then bound to change the nature of the many-body wave- 
function of the electron gas in the dot. For example, within the approximate Hartree-Fock picture, the effective 
single particle wave functions will be modified as electrons are kept added to the dot, resulting in the breakdown 
of the Koopmans picture ||l^. It was shown that the Koopmans picture, which asserts that the effective single- 
particle Hartree-Fock states remain unchanged as electrons are added to or removed from the system, breaks 
down when it concerns with finite (and disordered) systems, with sufficiently strong electron-electron interaction. 
For the systems studied that breakdown happened at r^ < 1.5 where rs = Y^3/(47rnoao^) with no the electron 
density and gq the Bohr radius. A remarkable experimental evidence has been provided by the experiment of 
[pO|, where the magnetofingerprints of various excited states of a dot at different particle number have been 
compared. It turned out that by adding electrons to the dot of concern (the total number of electrons was 
200) the magneto fingerprints have been significantly modified, indicating scrambling of the electronic states. 
A more systematic study of this scrambling has been reported in ref. ||l^ . We speculate that the suppression 
of the conductance correlations in the present context may be another tool to evaluate the scrambling. One 
systematic measurement which is called for is to repeat the experiment for different values of rs . Evidently, to 
facilitate a more detailed comparison with our calculated expression one would need data concerning various 
values of U/A, varying temperature and a larger number of valleys. 
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APPENDIX A: A FOR FINITE TEMPERATURE 



For the conductance correlation function we need the integral 



A{Ni,x,N2,x) = duj 



/(cat - w) 



ejv — tj + xU + iV cj — ejv + (1 — x)U — iV 



eN + d — LU + xU — iT lj — ejv — d + (1 — x)U + iV 



(Al) 



where d = sn^ — e^i ■ Since we are primarily interested in the correlation between different valleys we set xi = X2 = x. 
By means of contour integration we obtain 

A - rf(Y3^5^p^[V'(^ - Pi^U + d)) - 4'il - pxU) + ^(i - m - x)U + d)) - ^{^- - /3(1 - x)U)] + 
- id + u)ii - eM-m '^^^ - ^^' + "^^^ ' ^^^ + ^^' ' "^^^ + ^^^ " ^^' + ^' ' "^^^^ ' ^4 + ^"^^^ ' 



where -0 is the Digamma function and /? = j3/{2'Ki). For the normalization in (18) we need A for d -^ 0. We obtain 
A(iV, X, TV, x) = * (" A _ ^"j ("^(1) (i _ fiu{l -x))+ ^(1) (1 - /3C/.X)') - 
- ^ (^^'^(^ - /3C/(1 - x)) + 0(^)(i - /3C/X)) 
with 0(^'^) the first and second derivative of the Digamma function. 



(A3) 



APPENDIX B: A TOY MODEL FOR THE DISORDER 



We employ here a toy model which facilitates the calculation of the correlation function C yet, we feel, captures some 
of the essential physics. This model has been previously employed in ref. |^ as part of our analysis of the transmission 
phase correlations through a quantum dot. The model is based on the following simplifying assumptions: 

(i) The level spacing is constant, i.e. e^ = A • j. This mimics the level repulsion typical for a disordered system. 

(ii) For the product of tunneling matrix elements we assume V^^V^ = V^ai with a^ a random variable which can 
take the values +1 and —1 with equal probability. This models the fluctuations in the wave function caused by 
the disorder. 



With this model we calculate the correlation function of transmission amplitudes in the cotunneling regime 

Ct{N^,xi,N2,X2) = {t{Ni,xi)t{N2,X2r) (Bl) 

and show that we obtain the same results as with the formalism of section O. For one disorder configuration (a 
sequence of the Uj) the transmission amplitude in the cotunneling regime (with N electrons) at low temperature 
kT <^ A becomes 



'-"MEijT 



E 



E- {{k - M)A - xU) + iT 



]>M+ 



E-{{k- M)A + U{1 - x)) + il 



(B2) 



The disorder average is now equivalent to an average over the at. We find (a^) = and (aiaj) = 6ij. We set £^ = 
and obtain 



y2 /^ 1 

C.iN,x,,N + M,X2) = ^(Y.j-j^-^^ 



1 



Uk-N -M -X2U 



N+M 



■^ 1 1 

^ ^ k-N +U{l-xi)k- N - M -X2U 



E 



1 



1 



k - N +U{1 - xi) k - N - M +U{1 - X2) 



(B3) 



fe=Ar+j\/+l 

where hi — U/A. The sums can be expressed in terms of the Digamma function ^ 

M-l 



E 

n=0 



1 



1 



{n + a){n + b) 



a 



-iijib) - yj{a) + ^{M + a) - V(Af + b)) 



(B4) 



We arrive at 



^C.(iV,.,A. + M,..)^^^^^^-l--^ 



[tP{M+Ux2)-ij{Uxi) 



+i^{N + Uxi) - tP{N + M+UX2) + V^(l + M+U{1 - xi)) - V(l +Z^(1 - X2))] 

+ Tmr-^ ^n [^(l + ^^ + "(l - ^i) - ^-(1+^(1 - a^i) + V^(M + Ux2 - mx2)] 

M + (1 + 2:2 — xijU 



(B5) 



When we are sufficiently far away from the resonances xi^2hl ^ 1 we can use the asymptotic expansion of the Digamma 
function 4'{z) ~ logz for z ^ 1. In the limit iV — > oo we obtain 



^Ct{N,xuN + M,X2) = j^r-T^ ^ 

K^ M + [X2 — xi)U 



log 



M + X2Ul + M+{\-xi)U 



+ 



1 



M+{l+X2-xi)U 



log 



xiU 
M + x2Ul + M+{l-xi)U 



1 + (1 - X2)U 



X2U 



l + {l-xi)U 



With the formalism of section || the correlation function Ct is given by 



Ct{Ni,xuN2,X2) = A\g''%g'^A{Nuxi,N2, X2) 



(B6) 



(B7) 



with A defined in ( |22| ) and calculated in (|2j) . Comparison with ( pq ) shows that 
T/2 = A2(?igi)/e2(ag«)/e2. 
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